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I ABSTRACT 

We use photometric data from the Spitzer mission to explore the mid- and far-infrared prop- 
erties of 10 red QSOs (J-K>2, R-K>5 mag) selected by combining the 2MASS in 

I the near-infrared with the SDSS at optical wavelengths. Optical and/or near-infrared spectra 

are available for 8/10 sources. Modeling the Spectral Energy Distribution (SED) from UV to 
far-infrared shows that moderate dust reddening (Ay = 1.3 — 3.2) can explain the red op- 
tical and near-IR colours of the sources in the sample. There is also evidence that red QSOs 
have 60/12/im luminosity ratio higher than PG QSOs (97 per cent significance). This can be 

I interpreted as a higher level of star-formation in these systems (measured by the 60 /im lumi- 

nosity) for a given AGN power (approximated by the 12 luminosity). This is consistent 
with a picture where red QSOs represent an early phase of AGN evolution, when the super- 

I massive black hole is enshrouded in dust and gas clouds, which will eventually be blown out 

(possibly by AGN driven outflows) and the system will appear as typical optically luminous 
QSO. There is also tentative evidence significant at the 96% level that red 2MASS QSOs are 
more often associated with radio emission than optically selected SDSS QSOs. This may indi- 
cate outflows, also consistent with the young AGN interpretation. We also estimate the space 

I density of red QSOs relative to optically selected SDSS QSOs, taking into account the effect 

of dust extinction and the intrinsic luminosity of the sources. We estimate that the fraction 
of red QSOs in the overall population increases from 3% at Mk ~ —27.5 mag to 12% at 

I Mk — —29.5 mag. This suggests that reddened QSOs become more important at the bright 

end of the Luminosity Function. If red QSOs are transition objects on the way to becoming 
typical optically luminous QSOs, the low fractions above suggest that these systems spent less 
than 12% of their lifetime at the "reddened" stage. 
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1 INTRODUCTION 

In the last few years an increasing body of evidence points to an 
intimate link between the formation of galaxies and the growth of 
the supermassive black holes (SBH) at their centres. For example, 
the accretion history of the Universe (Barger et al. 2005; Hasinger 
et al. 2005) looks striki ngly sim ilar to th e evolution of the cosmi c 
star-formation rate (e.g. lHopkin s 2004; H opkins & Beacomll2006h . 
Additionally, most ne arby spheroids contain a SBH (SBH; e.g. 
iMagorrian et ai] |l998l) . with a mass that scales with the stellar ve- 
locity dis£ereion_oftheJiost galaxy bulge (e.g lFerrarese & MerrittI 
|200(]| ; lOebhardt et al.ll200Ci|) , suggesting similar formation epochs 
and/or mechanisms. 

Recent cosmological simulations account for the observations 



above by assuming that both the gal axy and the SBH form during 
galax y interactions and merge rs (e.g Kauffmann & Haehnelt 2OO0l 
2002i ; iDi Matteo et"ai] |2003| ; iMenci et al.. ,2004; Cattaneo et aU 



20061) . Such violent events are very efficient in driving gas to the 



nuclear regions, where it forms stars or it is consumed by the SBH. 
Because the nuclear regions of the interacting system are dusty, 
both processes above most likely take place within a cocoon of dust 
and gas clouds. Energetic considerations underline the significance 
of the energy released by th e AGN in the process above, which can 
regulate the gas inflows (e.g lSilk & Reeslll 9981 ; [pabianlll 9991 ; iKind 
|2003). Simulations of major mergers indeed show that the feedback 
from the AGN can be particularly violent close to the final merging, 
leading to the blow out of the dust and gas clouds and the termina- 
tion of both the star-formation and eventual l y the AGN activity (e .g. 
iDi Matteo erai]|2005l ; [Springel et alj2003 ; lHopkins et al.ll2006l) . 
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The key prediction of the picture above is that there 
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should be an association between SBH accretion, starbursts and 
mergers. Although there are examples of s uch systems (e.g. 



Genzel et al. 1998; Cid Fernandes et al. 2001; Francesc hini et al.l 
2003; Georgakakis et al. 2004; Alexander et al. 2005), they ap- 
pear to be the exception rather than the rule. X-ray surveys, ar- 
guably the most efficient method for finding AGN, have shown that 
the majority of X-ray s ources are hosted by red evolved galax- 
ies I^Na ndra et J] l2007h . with prominent bul ges and little mor- 
phological eviden ce for ongoing interactions tarogin et alj|2003l ; 
IPierce et al.ll2007h . This may suggest that X-ray wavelengths, be- 
cause of their high sensitivity to AGN, select a large number of 
systems after the peak of their activity (i.e. final merging) and dur- 
ing the decline stage of the SBH accretion. In order to understand 
the interplay between SBH and galaxy formation one should study 
the properties of AGN hosts (e.g. star-formation rate, optical mor- 
phology) close to the peak of the activity, when the system appears 
as luminous QSO. Such studies however, are hampered by the enor- 
mous brightness of the QSO nuclei relative to the host galaxy. Nev- 
ertheless progress has been made in the last few years. Mid-Infrared 
spectroscopy with Spitzer reveals PAH features, i.e. evidence for 
star-formation, in about 30% of optically selected luminous QSOs 
dSchweitzerlT al. 2006; Netzeretal. 2007). HST high resolution 
imaging and ground based observations with adaptive optics find 
that a fraction of the QSO population is hosted by either ellipti- 
cal galaxies with fine structure indicative of past interactions (80%; 
I 1' ' 

IBenn ert et al. 2008), or galaxies with obvious signs of disturbance 
(30%; Guyon et al. 2006). The evidence above suggests that at least 
some QSOs are formed in major mergers accompanied by starburst 
events. 

An important development in studies of the interplay be- 
tween AGN and galaxies has been the identification of QSOs with 
optical and/or near-infrared (near-IR) colo urs much redder than 
those of typical UV/ optically selected OSOs. lCutri et al.l ( l200lh and 
IWilkes et alj ( |2002|) for example, have found a previously unde- 
tected population of broad line AGN among sources w ith J — K > 
2 mag in the Two Micron All Sky Survey (2MASS). Iwhite et al.l 
l l2003 h studied an 7-band magnitude limited (J < 20.5 mag) sam- 
ple of radio sources and showed that the QSOs in that sample have 
redder B — R colours, on average , than those selected at bluer 
optical filters. iGlikman et all ( |2004 l2007f) demonstrated that se- 
lecting radio (1.4 GHz) sources in the region of the colour space 
J—K > 1.7 and R—K > 4.0 mag yields a 50 per cent success rate 
of discovering QSOs that are substantially redder than those found 
in optical sur veys. Similar results were obtained more recently by 
lUmitia et alj ( l2008a ) using a slightly modified colour wedge for 
selecting radio sources ( J — K > 1.3 and r — K > 5.0 mag). 
The optical/near-IR continua of the red QSO population identi- 
fied in the studies above are con sistent with moderate amounts of 
dust reddening {Ay « 1 5 - 3 IWhite et ar]|2003l ; iGUkman et al.l 
I2OO7I ; lUrrutia et ai]|2008ah. A n intrinsically flat UV/optical spec- 
trum (e.g. [Richards et al H003I) or synchr otron emission with a high 
turnover frequency I Whiting et alj|200ll) can also explain the red 
colours without invoking dust. It is unlikely however, that these 
alternative scenarios apply to the entire population of red QSOs. 
It is therefore proposed that this new population of dust reddened 
QSOs represent systems shortly before or during the blow out stage 
of AGN evolution, when the central source is still cocooned in 
dust and gas clouds. If true, reddened QSOs are intermediate be- 
tween dusty starbursts and UV/optically luminous QSOs. Clues on 
the evolutionary stage of these systems can be obtained by study- 
ing the properties of their host galaxies. Deep optical imaging for 
example, has shown that the morphology of the red QSO hosts 



are in most cases disturbed, suggesting a higher fraction of ongo- 
ing interactions compared to typical UV/optically luminous QSOs 
( iHutchings etai]|20C>3l , l200^ ; IUrrutia et alj|2008bl) . consistent with 
the young AGN picture. 

Another property of red QSOs that has not yet been addressed 
and is particularly relevant to their evolutionary stage is the level of 
star-formation rate in the host galaxy. In this paper we address this 
issue using Spitzer to constrain and to model the mid- and far-IR 
SEDs of a sample of red QSOs identified in the Two Micron All 
Sky Survey (2MASS). Under reasonable assumptions, the shape 
of the IR SED is a powerful diagnostic of the physical processes 
(star-formation vs AGN act ivity) responsible for heating the dust 
in extragalactic sources (e.g. lRowan-Robinson et al]|2008l) . 

The fraction of red QSOs in the overall population is also an 
important parameter. In the scenario where these systems represent 
a stage of the AGN phenomenon, their fraction provides clues on 
the timescale of this stage. Currently, there is large uncertainty in 
the fraction of red QSOs missing from optically sele cted samples, 
with values ranging from 1 5 to over 50 per cent (e.g.lWilkes et al.l 
2002! ; [Richards et al]|2003l ; IWhite et al.ll2003l ; iGlikman et alj2004 



120071) . The uncertainty is associated with the selection of the ap- 
propriate comparison sample of optically selected QSOs and with 
difficulties in accounting for dust extinction in the selection of red- 
dened QSOs. In this paper we address these issues by using the 
1/Vmax formalism to estimate the maximum volume that a source 
is detectable, given the survey limits and taking into account the 
effect of dust extinction and the intrinsic luminosity of the source. 
For comparison we use the Slo an Digital Sky Survey (SPSS) op- 
tically selected QSO sample of ISchneider et al.l ( l2005h , which is 
large and with well defined selection criteria. Throughout the paper 
we adopt Ho = 70kms~^ Mpc"\ = 0.3 and SIa = 0.7. 



2 SAMPLE SELECTION 

Candidates for dust enshrouded active galaxies were selected in the 
All Sky Data Release of the 2MASS Point Source Catalog (PSC) 
by cross-correlating the positions of the 2MASS sources with the 
3rd data release (DR3) of the SDSS. 

First, 2MASS-PSC sources with Ks < 14.5 mag and J - 
Ks > 1.5 mag were selected. The magnitude limit is close to 
the 99 pe r cent completen ess limit of the 2MASS-PSC {Ks = 
14.3 mag; ICutri et al .I2OO5I), while the colour cut is to avoid Galac- 



tic stars (e.g. Francis et alj2004l) . We exclude sources for which the 
A's-band photometric quality flag (PH_QUAL parameter) has val- 
ues "X" and "U", which signify bad pho tometry and an up per limit 
in the Ks-hwA magnitude respectively jCutri et al]|2005l) . 

The resulting 2MASS-PSC source list was cross-correlated 
with the SDSS using a matching radius of 3 arcsec, which is much 
larger than the sub-arcsec positional accuracy of the SDSS and the 
2MASS surveys. For the SDSS magnitude limit of r ~ 23 mag we 
expect about 2 per cent spurious alignment rate within this search 
radius. The SDSS r-band AB magnitudes were converted to Vega 
i?-band magnitudes using the colour transformations of Fukugita 
et al. (1996). This is to select candidates for dust enshrouded ac- 
tive systems using the colour cut [R — Ks)vega > 5 mag, i.e. 
similar to th at of Extre mely Red Objects (EROs; Elston, Rieke & 
Rieke 1988). |Pozzetti & Mannucci (2000) have used the R- K vs 
J — K colour diagram to discriminate between early-type and dusty 
EROs. This is shown in Figure[T] We focus on the dusty subregion 
of this figure. A total of 10 sources of our SDSS/2MASS sample 
lie in the dusty part of the colour-colour space of Figure [T] These 
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Figure 1. J - K against R - K colour plot for the 2MASS/SDSS sample 
with K < 14.5 mag . The sample of dusty active galaxy candidates is 
shown with filled circles. The dashed lines show the colour criteria used to 
select these sources. The horizontal line corresponds lo R — K = b mag 
and the diagonal line has been introduced by Pozzetti & Mannucci (2000) 
to separate dusty from early-type Extremely Red Objects. 



are presented in Table [T] and form the main sample of this paper. 
The offsets between the 2MASS and the SDSS source positions are 
<0.4 arcsec for all 10 selected sources and therefore the probability 
of spurious identifications is small, < 0.04 per cent, for the SDSS 
magnitude limit of r ~ 23 mag. We note that the colour selec- 
tion criteria adopted here are similar to those used by ('Wilkes et al.l 
[2002, B - Ks > J - K > 2 mag) and ( Glikman et al. 2004. 
R-Ks> 4.0, J -K > 1.7 mag) to select 2MASS red AGN. 



3 OBSERVATIONS 

In addition to the available optical (SDSS; ugriz) and near-IR 
(2MASS; JHK) broad band data, which are presented in Table[T] 
we have also carried out mid- and far-IR photometry with Spitzer 
as well as optical and near-IR spectroscopy to determine redshifts. 
These observations are described below. 



3.1 Mid- and far-Infrared photometry 

The Spitzer mid- and far-infrared photometric observations were 
performed during Cycle-2 as part of the General Observer program 
with identification number 30607 (PI Georgakakis). The mid- and 
far-IR flux densities of the target sources are presented in Table [2] 
The 3.6 — 8.0/im data were taken with the IRAC (Infrared 
Array Camera) between July 9 and December 28 2006. Each ob- 
servation consists of 5 separate integrations of 2 s each organised in 
a random Gaussian dither pattern. The median separation between 
successive offsets is 0.9 arcmin. The observations were run through 
the S 14.4.0 version of the Spitzer Science Center pipeline described 
in IRAC data handbook. The basic calibration data (BCD) prod- 
ucts are processed for dark current, bias offset, linearisation, flat- 
fielding, cosmic ray detection and the flux calibration. The individ- 
ual frames of a given source were then registered and co-added into 



the final mosaiced image. The uncertainty in the calibration fac- 
tor applied to the final mosaics is about 10 per cent. Following the 
IRAC data handbook flux densities are estimated by using apertures 
of 6 arcsec radius and then applying aperture correction factors to 
account for the extended Spitzer point spread function (PSF) in dif- 
ferent wavebands. These are estimated 1.06 for IRAC 3.6, 4.5 /im 
bands and 1.10 for the 5.8 and 8.0 /im bands. Colour corrections 
are expected to be small as the mid-IR spectra of the sources in 
the sample approximate a power-law with ffv ~ constant. We 
therefore choose not to apply colour corrections. 

The 24, 70, and 160/im data were taken with MIPS (Multi- 
band Imaging Photometer for Spitzer) between July 14 2006 and 
September 15 2007. The observations used the MIPS small-field 
photometry mode with an exposure time of 3 s for individual 
frames. The 24 and 70/^m data were obtained in one cycle, while 
the 160/^m observations were carried out in 3 cycles. The total on- 
source exposure times were 48, 38 and 25 s at 24, 70, and 160/.im 
respectively. The MIPS images were created from ra w data frames 
using the MIPS Data Analysis Tools (MIPS DAT; 1 Gordon et all 
I2OO5I) version 3.10 along with additional processing steps. The pro- 
cessing steps for the 70 and 160/^m data are similar, but the steps 
for the 24 /^m data are significantly different from these other two 
bands. The 24 /im data processing is described first followed by 
descriptions of the 70 and 160/im data processing. 

The individual 24 /im frames were first processed through a 
droop correction (removing an excess signal in each pixel that is 
proportional to the signal in the entire array) and were corrected for 
non-linearity in the ramps. The dark current was then subtracted. 
Next, scan-mirror-position dependent flats were applied to the data, 
latent images were removed, and scan-mirror-position independent 
flats were applied to the data. Following this, the zodiacal light was 
measured by fitting planes to the background regions in each frame, 
and these planes were then subtracted from the data. Next, a robust 
statistical analysis was applied to cospatial pixels from different 
frames in which statistical outliers (e.g. pixels affected by cosmic 
rays) were masked out. After this, final mosaics were made with 
pixel sizes of 1.5", residual backgrounds in the images were sub- 
tracted, and the data were calibrated into astro nomical units. The 
calibration factor for the 24 /im data is given by Engelbra cht et aEI 
('2007) as (4.54 ± 0.18) x 10"^ MJy sr"^ [MIPS instrumental 
unit]^^. A 10.5 arcsec radius aperture is used for the photometry. 
The correction factor that accounts for PSF losses is estimated to be 
1.22. The mid-IR spectra of the sample sources is close to a power- 
law. For this SED type the colour correction at the 24 /im band is 
small (<3 per cent) and we choose not to apply it. 

In the 70 and 160/xm data processing, the first step was to fit 
ramps to the reads to derive slopes, during which readout jumps 
and cosmic ray hits were also removed and an electronic nonlin- 
earity correction was applied. Next, the stim flash frames taken by 
the instrument were used as responsivity corrections. The dark cur- 
rent was subtracted from the data, an illumination correction was 
applied, and short term variations in the the signal (often referred 
to as drift) were removed. Next, a robust statistical analysis was 
applied to cospatial pixels from different frames in which statisti- 
cal outliers (e.g. pixels affected by cosmic rays) were masked out. 
Once this was done, final mosaics were made using square pixels 
of 4.5 arcsec for the 70/im data and 9 arcsec for the 160 /im data. 
The residual backgrounds were measured and subtracted from the 
images. Finally, flux calibration factors were applied t o the data. 
The 70/im calibration factors given bv lGordonetd] 1200% are 
702 ± 35 MJy sr"^ [MIPS i nstrumental unit]~\ aii d the 160/im 
calibration factor is given bv lStansberrv et al] j2007h as 41.7 ± 5 
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Table 1. The 2MASS selected active dusty galaxy candidates. The columns are: (1): source name; (2): right ascension of the optical source; (3): declination of 
the optical source; (4): M-band AB magnitude and error; (5): g-band AB magnitude and error; (6): r-band AB magnitude and error; (7): i-band AB magnitude 
and error; (8): z-band AB magnitude and error; (9): J-band Vega magnitude and eiTor; (10); H-band Vega magnitude and error; (11): Ks-band Vega magnitude 
and error; 



ID 


a 


<5 


2MASS.02 


20h56ni29.76s 


-06d50m554s 


2MASS.03 


20h48ra37.25s 


-00d24m37.3s 


2MASS.04 


17h04ra2l.l6s 


+24d33m40.9s 


2MASS.05 


15h33m34.18s 


+07d35m05.7s 


2MASS.06 


13h40ra39.68s 


+05dl4ml9.9s 


2MASS.07 


13h21m22.89s 


+45d02m24.8s 


2MASS.08 


12h52ml2.93s 


+07dl5m04.7s 


2MASS.10 


08h25m2.06s 


+47dl6m52.0s 


2MASS.1 1 


01h34m35.67s 


-09d31m03.0s 


2MASS.12 


00h36ra59.82s 


-01dl3m324s 



H 



K 



21.75 
20.48 
22.40 
22.12 
20.69 
24.52 
21.92 
21.38 
25.52 
22.05 



± 0.36 
± 0.13 
± 0.42 
± 0.94 
± 0.10 
± 2.26 
± 0.33 
± 0.29 
± 0.81 
± 0.31 



20.34 
20.12 
21.57 
22.66 
19.82 
22.44 
21.30 
20.88 
23.51 
21.52 



± 0.04 
± 0.03 
± 0.07 
± 0.50 
± 0.02 
± 0.22 
± 0.06 
± 0.07 
± 0.32 
± 0.07 



19.59 
19.41 
20.23 
20.97 
18.97 
20.75 
21.03 
20.50 
21.20 
20.35 



± 0.02 
± 0.02 
± 0.03 
± 0.16 
± 0.01 
± 0.06 
± 0.07 
± 0.07 
± 0.05 
± 0.03 



18.69 
18.93 
19.36 
19.65 
18.46 
19.36 
20.30 
19.68 
19.60 
19.88 



± 0.01 
± 0.02 
± 0.02 
± 0.07 
± 0.01 
± 0.03 
± 0.05 
± 0.05 
± 0.02 
± 0.03 



18.47 ± 
18.32 ± 

18.65 ± 

18.60 ± 
17.70 ± 

18.61 ± 
19.43 ± 
19.56 ± 
18.31 ± 

18.66 ± 



0.04 
0.04 
0.05 
0.11 
0.02 
0.05 
0.09 
0.16 
0.03 
0.04 



< 16.48 
16.67 ± 0.15 
16.82 ± 0.17 
16.98 ± 0.23 
16.16 ± 0.11 
16.58 ± 0.13 

< 17.88 
17.14 ± 0.22 
16.18 ± 0.13 
16.56 ± 0.10 



15.54 
15.08 
15.72 
15.82 
14.80 
15.12 
16.53 
15.84 
14.79 
15.10 



± 0.09 
± 0.08 
± 0.13 
± 0.14 
± 0.08 
± 0.08 
± 0.18 
± 0.18 
± 0.07 
± 0.06 



13.76 
13.64 
14.49 
14.44 
13.25 
14.38 
14.49 
14.16 
13.58 
13.63 



± 0.04 
± 0.04 
± 0.08 
± 0.11 
± 0.03 
± 0.09 
± 0.10 
± 0.06 
± 0.05 
± 0.04 



Table 2. The 2MASS selected active dusty galaxy candidates. The columns are: 1: source name; 2: log of IRAC 3.6^m flux density in /iJy. The error includes 
both background and callibration uncertainties; 3; log of IRAC 4.5/jm flux density in /.tJy; 4: log of IRAC 5.8^m flux density in ^Jy; 5: log of IRAC S.O^m 
flux density in fiJy; 6; log of MIPS 2Afim flux density in /.iJy; 7: log of MIPS 70/im flux density in /.tJy; 8: log of MIPS IGOfim flux density in /.tJy; 9: 
1.4GHz radio flux density in mJy. All measurement are from FIRST except for sources 2MASS_2 and 2MASS_3, which are from the NVSS. 10: Redshift. 
Photometric redshifts are marked with "ph". See text for details. 



ID 


3.6 ^Ta 


4.5 


5.8 


8.0 Mm 


24 ^ra 


70 


160 


1.4 GHz 


red.shift 


2MASS.02 


3.91 ± 0.01 


4.10 ± 0.01 


4.29 ± 0.01 


4.40 ± 0.01 


4.52 ± 0.01 


< 4.57 


< 4.92 


5.6 


0.635 


2MASS.03 


3.93 ± 0.01 


4.04 ± 0.01 


4.19 ± 0.01 


4.28 ± 0.01 


4.35 ± 0.02 


5.71 ± 0.03 


5.32 ± 0.05 


3.6 


0.433 


2MASS.04 


3.60 ± 0.01 


3.74 ± 0.01 


3.85 ± 0.01 


4.01 ± 0.01 


4.26 ± 0.02 


5.32 ± 0.03 


5.32 ± 0.05 




0.422 


2MASS.05 


3.27 ± 0.01 


3.45 ± 0.01 


3.74 ± 0.01 


3.99 ± 0.01 


4.36 ± 0.01 


5.28 ± 0.03 


5.12 ± 0.05 




3.10^*^ 


2MASS.06 


3.96 ± 0.01 


4.09 ± 0.01 


4.25 ± 0.01 


4.36 ± 0.01 


4.78 ± 0.01 


5.41 ± 0.03 


5.09 ± 0.05 


8.9 


0.264 


2MASS.07 


3.40 ± 0.01 


3.62 ± 0.01 


3.83 ± 0.01 


4.03 ± 0.01 


4.29 ± 0.02 


4.77 ± 0.03 


< 4.51 




2 73Pfi 


2MASS.08 


3.08 ± 0.01 


3.26 ± 0.01 


3.44 ± 0.02 


3.74 ± 0.01 


4.39 ± 0.01 


5.26 ± 0.03 


< 4.85 




2.160 


2MASS.10 


3.80 ± 0.01 


3.98 ± 0.01 


4.14 ± 0.01 


4.28 ± 0.01 


4.54 ± 0.01 


4.96 ± 0.03 


< 4.95 


63.2 


0.804 


2MASS.1 1 


3.53 ± 0.01 


3.67 ± 0.01 


3.91 ± 0.01 


4.07 ± 0.01 


4.18 ± 0.02 


< 4.57 


4.51 ± 0.05 


919.7 


2.210 


2MASS.12 


3.89 ± 0.01 


4.02 ± 0.01 


4.17 ± 0.01 


4.29 ± 0.01 


4.78 ± 0.01 


5.82 ± 0.03 


5.44 ± 0.05 


0.8 


0.291 



MJy sr~^ [MIPS instrumental unit]~^. For the photometry we use 
apertures with radii 4.5 and 9 arcsec at 70 and 160fim respectively. 
In this case, the correction factors that account for PSF losses are 
1.30 and 1.87 respectively. Colour corrections for the MIPS 70 and 
IGOfim bands are typically smaller than the calibration uncertain- 
ties and are therefore not applied to the data. In the case of non- 
detection we assign an upper limit to the flux density which corre- 
sponds to 5 times the standard deviation of the background. 

3.2 Spectroscopy 

Optical spectroscopy for the target sources was carried out at the 
Kitt Peak National Observatory (KPNO) 4-m telescope and the 
4.2m William Herschel Telescope (WHT). 

The KPNO observations used the Ritchey-Chretien spectro- 
graph in single-slit mode with the BL 181 grism blazed at 7500A. 
This setup provides a resolution of about 7A in the wavelength 
range 5500 — 9500A. An exposure time of 20min was adopted. 
These observations were carried out in March 10 2006. 

The WHT spectroscopy used the ISIS (Intermediate disper- 
sion Spectrograph and Imaging System) double armed spectro- 
graph during service time in June 17 and July 8 2007. The obser- 
vations were made using the 5300A dichroic and the R300B and 
R158R gratings in the blue and red arms of the spectrograph re- 
spectively. The spectral resolution was about 7A for the red arm 
and 4A for the blue arm spectrum. The total on-source integration 
time was 20min split into two lOmin exposures. 

The data were reduced following standard methods as imple- 
mented in the NOAO package of IRAF. After subtraction of the bias, 
the data were flat-fielded using observations of internal calibration 



lamps. Cosmic ray events were removed usin g the Laplacian Cos - 
mic Ray Identification package (LACOSMIC;|van Dokkunjl200lh . 
The one dimensional spectra were extracted and wavelength cal- 
ibrated using observations of CuNe and CuAr arc lamps. For the 
flux calibration observations of the spectrophotometric standard 
stars BD 282411 and BD 332642 were used. Because of the ap- 
parent faintness of our targets at wavelengths shorter than 5000A 
the WHT/ISIS blue arm spectra did not show any signal and are 
therefore not presented in this paper. 

Near-infrared spectroscopy of SDSS/2MASS EROs was ob- 
tained in queue mode on the UK Infrared Telescope (UKIRT). The 
observations used the UKIRT Imaging SpecTrometer (UIST) in 
long slit mode with a mixture of the IJ and HK grisms, the exact 
selection depending on the expected redshift of the object based 
on photometric redshift estimation. The observations were carried 
out in clear conditions on 18, 22, 23 June, 2007. Integration times 
were ~ 1.5 hours for the HK grism and '-^ 45 minutes for the IJ 
grism. The data were reduced using the automated ORACDR sys- 
tem which performed flat fielding, wavelength calibration based on 
arc spectra, subtracted the chopped images and coadded separate 
sub-exposures. Flux calibration was carried out using the BS5019 
and BS8 spectrophotometric standards. No correction for slit losses 
were applied so the flux calibration is correct in a relative sense. 
The spectra were then analyzed using a combination of IRAF and 
IDL routines to determine redshifts and to extract relevant emission 
line parameters. 

Spectroscopic redshifts are available for 8/10 2MASS sources, 
7 from our own observations and 1 from the literature. The lack of 
spectroscopic redshift determinations for 2 sources in the sample is 
because of low S/N optical spectra. The optical and near-IR spectra 
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of the sources with successful redshift determinations are shown 
in Figure |2l The redshifts are presented in Table |2] and estimates 
of the FWHM (Full Width Half Maximum) of emission lines are 
discussed in the Appendix and are listed in Table[3]. 



4 SPECTRAL ENERGY DISTRIBUTIONS 

The observed optical to mid-IR Spectral Energy Distribution (SED) 
of the sampl e sources are modeled f ollowing the methods fully 
described in iRowan-Robinson et all ( |2005|, |2008|) . In brief the 
[/-band to 4.5 /im photo metric data are fit us ing a library of 
8 templates described by iBabbedgeetd] ( l2004h . 6 galaxies (E, 
Sab, Sbc, Scd, Sdm and sb) and 2 AGN. At longer wavelengths 
(5.8 — 24 /im) any dust may significantly contribute or even dom- 
inate the observed emission. Before fitting models to these wave- 
lengths the stellar contribution is subtracted from the photomet- 
ric data by extrapolating the best-fit galaxy template from the 
previous step. The residuals are then fit with a mixt ure of four 
templates: cirrus Jifstathiou & Rowan-Robinso r]|2003h. AGN dus t 
tori jRowan-Robiriso n 1995; Efstathiou & Rowan-R obinsonI 1995h . 
M82 and Arp220 starbursts (Efstathiou et al. 2Q0(Q. The model- 
ing above provides both information on the dominant emission 
mechanism in the optical and the infrared (AGN vs star-formation) 
and an estimate of the total infrared luminosity, Ltot, of our 
sample so urces in the wavelengt h range 3 — lOOO^m. As dis- 
cussed by iRowan-Robinson et all ( FOOS.) Ltot is expected to be 
accurate within a factor of two. In this exercise the redshift is 
fixed to the spectroscopically determined value, if available (8/10 
sources). For the two sources in the sample without spectroscopic 
redshifts we also estimate photometric redshifts, although we cau- 
tion that these are likely to be uncertain because of their extreme 
optical/near-IR colours. From the spectroscopic sample we esti- 
mate that the fraction of catastrophic redshifts, defined as those 
with {zspcc - Zphot)/{l + Zspec) > 0.15, is 38 per cent (3/8). 
The rms value of the quantity {zspec — 2phot)/(l + 2spec), after 
excluding catastrophic failures, is 0.07 and provides an estimate of 
the accuracy of the photometric redshifts. 

Figure [3] overplots the best-fit models to the observed SEDs 
of the red 2MASS sources. The derived parameters are presented 
in Table [3] In summary, the mid-IR SED of all sources is domi- 
nated by hot dust, which is fit by an AGN torus component. An 
additional starburst template is required to fit the far-IR data of 
some sources in Table [S] The optical part of the SED is fit with 
a reddened QSO tem plate adopting the SMC extinction curve of 
[Richards et sl\ j2003h . The derived optical extinctions are in the 
range Ay = 1.3 - 3.2 {E{B - V) x 0.4 - 1.1). 



5 NATURE OF THE 2MASS SELECTED SAMPLE 

The nature of the red 2MASS sources is explored by combining 
information from the optical/near-IR spectroscopy and the broad- 
band UV-to-far-IR photometry. There is strong evidence that all 
sources in the sample are powered by accretion on a central super- 
massive black hole. Firstly, for seven out of the eight sources for 
which optical and/or near-IR spectra are available, either from our 
own observations or from the literature, broad Balmer or Paschen 
emission lines are found with widths > 1000 km/s. For more de- 
tails on the spectroscopic properties of individual objects see the 
Appendix section. Additionally, the red 2MASS sources are lu- 
minous in the mid-IR and their SEDs at these wavelengths are 



fit by the hot dust AGN tori mode l s of Rowan-Robinsonl ( Il995h 
and lEfstathiou & Rowan-Robinsonl ( Il995[) . The total IR luminos- 
ity (3 — lOOO/im) of this component is estimated in the range 
IO12 _ iQi^Lg (4 X lO"*^ - 4 X 10*'' erg s"^), i.e. exceeding 
the limit for either ULIRGs (Ultra-Luminous Infrared Galaxies) or 
HyLIRGs (Hyper-Luminous Infrared Galaxies). A dopting an AGN 
bolometric correction factor of Lboi /Lm ~ 3 (e.g lRisaliti & ElvisI 
'2004*), the luminosity interval above exceeds the limit Lboi }t 
10'"'erg/s, which is often used to differentiate luminous QSOs 
from Seyfert galaxies. Therefore, in the following we refer to the 
sample of sources in Table[T]as 2MASS or red QSOs. 

The red optical/near-IR colours of the 2MASS QSOs in Table 
1 ca n be attributed to eith er an intrinsically red optical continuum 
(e.g. [Richards et al]|2003h . synchrotron emission w ith a turnover 
frequ ency in the UV/optical part of the spectrum (e.g. lWhiting et al.1 
I2OOII) or dust reddening. Each of these possibilities is discussed 
below. 

Radio selec ted QSOs from the Parkes half-Jansky flat- 
spectrum sample jPrinkwater et al.l 19971) have a large spread in the 
op tical/near-IR colours , with the redde s t obje cts having B — K > 
7 ( lFrancisetal1l2000l) . IWhiting et all (1200 1*) proposed that syn- 
chrotron emission with a turnover frequency in the range 0. 1 — 3/im 
can explain the red colours of a sizeable fraction of this population 
(about 40 per cent). This is relevant to this study as many 2MASS 
QSOs (6/10; see Table 1) are detected in the NVSS dCondon et al.l 
Il998i limit ^ 2.5 mjy) or the FIRST ("Becke r et alj|l995l, limit 
« 1 miy) radio ( 1 .4 GHz) surveys with luminosities in the range 
~ 10'^* — 10^'' W/Hz. Therefore, a contribution by synchrotron 
emission in the optical/near-IR wavebands is a possibility. Figure|4] 
explores this scenario by extrapolating the observed 1.4 GHz flux 
density to the UV assumin g a standard s ynchrotron spectrum of 
the form f^ocv "° (e.g. lCondon|[T993) . This figure shows that 
the observed radio emission is too faint to explain the observed 
red colours of most radio detected 2MASS QSOs, in agreement 
with previous studies (e.g. Glikman et al. 2007). A possible ex- 
ception i£_source_2MASS^l 1. This lensed source was first identi- 
fied by iGregg et alj ( l2002h because of its unusually red B ~ K 
colour. The radio flux density of this system is sufficiently bright 
(Ri 1 Jy) that can potentially affect the optical colours as suggested 
by Siting et al. (2 001J. 

Richards et al. showed the observed optical colours of 



QSOs are a strong function of redshift and proposed the use of 
relative colours (e.g. A(g — i); difference between the observed 
and the median colour of QSOs at the same redshift) to discrim- 
inate between intrinsically red and dusty systems. These authors 
showed that the A(g — i) distribution of SDSS optically selected 
QSOs can be approximated by a Gaussian, suggesting a range of 
intrinsic optical continuum slopes, and a red tail, which is domi- 
nated by dust reddened systems. All the 2MASS QSOs presented 
in this paper have l\{g — i) > 0.8, which places them in the red tail 
of the SDSS optically se lected QSO distribution (e.g. see Figure 6 
of [Richards et ai]|2003l) and in the region of the parameter space 
where sourc es with dust reddenin g E[B — V^) -SJ 0.1 are expected 
to be found. [Richards et alj ( 12003 ) shows that such red colours are 
most likely because of dust rather than a steep optical power-law 
spectrum. 

Dust reddening indeed provides good fits to the UV-to-near- 
IR SEDs of the 2MASS QS Os, if the SMC extinction curve is 
adopted dRichards et"ai]|2003l) . This type of extinction is strongly 
supported by the shapes of the observed UV/optical SEDs, which 
are smooth and do not show the characteristi c bump at rest-fram e 
fa 22OOA of the Galactic extinction curve ^Gordon et ai]|2003h . 
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Figure 2. Optical and near-IR spectra of the sources with successful redshift determination from our own spectroscopic observations. TTie gap between near-IR 
spectra corresponds the gap between the H and K windows of the HK grism. 
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Figure 3. Observed Spectral Energy Distiibutions of the red 2MASS sources. The dots are the observed UV-to-far-IR photometry and the continuous hnes 
are the best-fit models. Different sources are labeled and are offset by an arbitrary constant for clarity. The mid-IR wavelengths of all sources are fit by a hot 
dust torus component (dashed line) which dominates at these part of the SED. Some sources also require an additional starburst component (dotted line at 
infrared wavelengths) to fit the far-IR data. The optical part of the SED is fit with an reddened QSO template. The reddening is listed in Table[3] Some sources 
require an additional component to fit the u and/or g-band SDSS photometry (2MASS_3, 2MASS_10, 2MASS_12). A star-formation component is used here, 
although we cannot exclude the possibility of scattered AGN light. 



The derived optical extinctions for the 2MASS QSOs are in the 
range Ay = 1.3 — 3.2 (see TableO, much higher than optically se- 
lected QSO samples, which are typically sensitive to systems with 
Ay < 0.8 or E(B - V) < 0.2 for SMC type extinction (e.g. 
[Richards et aiT2003) . 

An independent estimate of the reddening can be obtained us- 
ing the Balmer or Paschen line decrements. For the former we as- 
sume Case-B recombination ( Storev & Hummer 1995), while for 
the latter we use empirically determined ratios, as Case B recombi- 
nati on is not a good appr oximation for the Paschen lines of QSOs 
(e.g lGlikman et al. 2007: Soiferetai]|2004h . This excercise is pos- 
sible for 4 out of 10 2MASS QSOs (see Appendix). These val- 



ues typically do not agree with the reddening derived from SED 
fitt ing. A similar trend using a larger sample has been reported 
bv lGlikman et al.l ( l2007h . This implies that the dust distribution is 
complex and therefore the line emission and the continuum are ab- 
sorbed by different material. 

The ratio of the total infrared luminosity in the dust torus 
(3 — lOOO^m) to the reddening corrected optical (0.1 — 3/im) 
QSO luminosity, Ltor/ Lopt, has been interpreted as a measure of 
the dust covering factor in AGN (Rowan-Robinson et al. 2008). For 
optically selected QSOs in the SWIRE survey the distribution of the 
\og{Ltor / Lopt) is well fitted by a Gaussian with mean -0.10 and 
standard deviation 0.26 (Rowan-Robinson et al. 2008). We find that 
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Table 3. The rest-frame properties of the 2MASS QSOs. The columns are: 1: name; 2: total infrared luminosity in solar units (Lq = 3.83 X 10^^ ergs~^) 
and in the wavelength interval 3 — 1000/^m. This is estimated from the SED fits described in section|4l 3: Infrared templates used to fit the mid- and far-IR 
observations. There are sources where both an AGN torus and a starburst component are required. 4: Ratio of the total infrared luminosity in the dust torus 
(3 — 1000 fim) to the reddening corrected optical (0.1 — Sfim) QSO luminosity. This ratio has been interpreted as a measure of the dust covering factor 
in AGN (Rowan-Robinson et al. 2008). 5: Optical V-band extinction in magnitudes estimated from fitting the observed SED as described in section |4] 6: 
Ks-bmid absolute magnitude. The Ks-hand flux density is estimated directly from the observations by simply convolving the rest-frame SED of each source 
(i.e. corrected for redshift) with the Ks-band filter function. 7: Radio power at 1.4 GHz in units of W/Hz. A radio spectrum of the form oc was 
adopted for the calculation of the k-correction. 8: FWHM in km/s of the broad component of the observed emission lines. We do not list FWHM for sources 
with narrow emission lines only (e.g. 2MASS_04) or without optical spectroscopy from our own follow-up program. The line used to measure the FWHM is 
also listed in brackets. 
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6/8 2MASS QSOs with spectrocopic redshifts (see Table 3), have 
\og{Ltor / Lopt) > 0.0, suggesting covering fractions larger than 
the typical SWIRE optically selected QSO. In this comparison we 
have excluded the two 2MASS QSOs without spectroscopic red- 
shifts to avoid uncertainties in the determination of their rest-frame 
properties because of errors in the photometric redshift estimates. 

In addition to the AGN activity, a number of sources in the 
sample also show evidence for star-formation. There are objects 
(see Table [3]l with excess emission in the far-IR over the ex trap- 
olation of the dust torus mode ls of [ Rowan- Robinsod h995l) and 
[sfstathiou & Rowan-Robinson In our modeling of the ob- 

served SEDs this excess is fit with cool dust associated with star- 
bursts. The total IR luminosity of this component is large, in the 
range ~ 10^^ — 10^^ Lq suggesting high star-formation rates, 
> lOOM0yr~^. Additionally, there are sources in the sample 
with excess emission in the bluer SDSS optical bands, g and/or 
u, over the expectation of the reddened QSO template (2MASS_3, 
2MASS_10, 2MASS.12). This excess can be fit by an additional 
starburst component that dominates over the reddened QSO emis- 
sion in the UV part of the spectrum. We cannot exclude however, 
the possibility of scattered radiation from the AGN itself as the ori- 
gin of this excess. It is interesting nevertheless, that 2MASS_3, one 
of the objects that shows the UV excess, has narrow emission lines 
in addition to the broad wings of the Ha. The flux ratios of the 
narrow emission line components place this source in the region 
of Tra nsition Objects in the diagnostic diagram of iKewlev et al.l 
suggesting some level of star-formation. 

In summary, the red 2MASS sources studied in this paper are 



dust reddened QSOs, while a large fraction (7/10 in Table[3l( have 
far-IR properties consistent with star-formation in the host galaxy. 

6 COMPARISON WITH OPTICALLY SELECTED QSOS 
6.1 Star-formation rate 

In this section we first explore possible systematic differences in 
the IR SEDs of the red 2MASS QSOs with optically selected ones. 
For this comparison we use the Palomar-Green (PG) QSO sample 
ISchmidt & Greei] ( Il983h . for which a rich set of multiwavelength 
observations is available. We consider in particular the sub-sample 
of 64 PG QSOs with mid- and far-IR observa tions (5 - 200/jm ) 
in the ISO archive analysed and presented by iHaas et al] ( l2003h . 
From that sub-sample only sources with SDSS optical photometry 
(total of 26) are used here. The near-IR data for these sources are 
from the 2MASS survey. In order to avoid uncertainties associated 
with the photometric redshift estimates, in this section we consider 
2MASS QSOs with spectroscopic redshifts only. The median red- 
shift of the PG QSO subsample used here is ~ 0.161, lower than 
that of 2MASS QSOs, z ^ 0.5. 

Figure [5] compares the rest-frame UV-to-IR SEDs of the 
2MASS QSOs with those of the optically selected PG QSOs. The 
SEDs are normalised to the rest-frame 12/im flux density, as this 
wavelength is proposed to be nearly u nbiased census of the AGN 
power (e.g. ISpinoglio & MalkaiJll989j ). As expected dust redden- 
ing makes the UV/optical SEDs of 2MASS QSOs steeper than 
PG QSOs. Also, there is evidence that the far-IR properties of 
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Figure 4. Spectral Energy distributions of the 2MASS QSOs with radio counterparts in the FIRST or the NVSS. The dotted Hnes coiTesp ond to the radi o flux 
density extrapolated to optical wavelengths assuming a radio spectrum of the form oc u-O-'^^ i.e. typical of synchrotron emission jCondorj[T993) . The 
radio flux densities are too faint to affect the optical/near-IR colours of the 2MASS QSOs. The o nly exception is 2MA SS_1 1 with a radio flux density of about 
1 Jy. Synchrotron emission with a turnover frequency in the UV/optical part of the spectrum (e.g. lwhiting et alj200ll) could be responsible for the red colours 
of this source. We note however, that the possibility of dust red dening cannot be ruled out. For ex ample, previous studies on the nature of this source have 
argued for dust associated with either the source or die lens fe.g. lHall et al.l2002l ; lGregg et alj2003) . 



the two samples differ. The normalised rest-frame far-IR flux den- 
sity of at least 4 out of 8 2MASS sources with spectroscopic red- 
shifts is elevated compared to the median SED of PG QSOs. Our 
SED modeling (s ee Table [3 i and r ecent studies of optica lly se- 
lected QSOs fe.g iHaas et aLlbOOl [Schweitzer et alJl2006h . sug- 
gest that the far-IR emission of these monsters is dominated by 
star-formation. Under this assumption, Figure |5] therefore shows 
that at least some 2MASS QSOs have, on average, higher star- 
formation than optically selected QSOs. This is further demon- 



strated in Figure |6] which plots 60 /im luminosity against 12 /im 
luminosity (see figure caption for details on the determination of 
the luminosities). For given 12 fim luminosity, some 2MASS QSOs 
are more luminous at 60 /im. The assessment of the significance 
of this excess is limited by the small sample size. Nevertheless, 
we use survival statis t ics as implemented in the ASURV pack- 
age dlsobe et al.lll986l : iLavallev et alJll992h to estimate the null 
hypothesis probability that the distribution of the luminosity ratio 
log i^L^(12/im)/!/L^(60/im) of 2MASS and PG QSOs is drawn 
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Figure 5. The rest-frame SEDs of 2MASS QSOs ai'e shown with the black 
dots connected with lines. Only sources with spectroscopic redshift are 
plotted. The long-dashed thick red line is the median rest-frame SED of 
optically selected QSOs. The red dotted lines correspond to the 25th and 
75th percentiles of the distribution. Survival statistics as implemented in 
the ASURV package (Isobe et al. 1986; Lavalley et al. 1992) are used to 
estimate the median and the 25th/75th percentiles of the PG QSO SEDs 
taking into account upper limits in the flux density measurements of in- 
dividual sources. All SEDs are normalized to the rest-frame 12/xm flux 
density, which provides a good proxy to total AGN power independent of 
dust extinction. 2MASS QSOs have much steeper SEDs in the UV/optical, 
as expected from dust reddening. Also some 2MASS QSOs have higher 
far-IR flux density (normalised to 12/xm flux density) compared to the me- 
dian SED of PG QSOs. The 2MASS QSO 2MASS.1 1 has significantly 
different SED from the rest of the 2MASS QSOs in the sample and is 
marked wi th the short-dashed black line. This is a lensed radio source at 
z RJ 2.2 jGregg et al.ll2002l) . Synchrotron em ission may affect th e rest- 
frame UV/optical SED of this radio source (e.g. lWhiting etalj200lh . 

from the same parent population. Using the Gehan's generalised 
Wilcoxon test we estimate this probabihty to be 3 per cent. The 
null hypothesis can therefore be rejected at the 97 per cent level 

(« 2a in the case of Gauss d istrib ution). 

ISchweitzer et alj ( l2006h and iNetzer et all jlOOTh find that a 
substantial fraction of optically selected PG QSOs show evidence 
for star-formation (up to 70% based on the far-IR detection rate, 
30% based on PAH detection rate). Moreover, they found that indi- 
cators of star-formation rate (e.g. PAH or 60/im luminosities) cor- 
relate well with indicators of AGN power (e.g. 6/im or 5100 A lu- 
minosities) for these systems. The red 2MASS QSOs studied here 
also show evidence for ongoing star-formation: 6 out of 8 sources 
(75%) with spectroscopic redshifts are fit by an additional starburst 
component in the IR (see Table O. An additional interesting re- 
sult from our analysis is that at the 97 per cent significance level 
the 2MASS QSOs have excess star-formation for their AGN power 
compared to the optically selected PG QSOs. 

6.2 Radio detection rate 

An interesting trend in Table 1 is the high fraction of the 2MASS 
QSOs with detection at 1.4 GHz: 6/10 are radio sources in either 
the NVSS of the the FIRST surveys. This is unexpected because the 



Figure 6. 60/im luminosity plotted against 12/im luminosity. Only 
2MASS QSOs with spectroscopic redshift are shown. The rest-frame 60/im 
and 12/im flux density is estimated from the observed SED by linearly 
interpolating between data points. Two out of eight 2MASS sources with 
optical spectroscopy (2MASS.08 and 2MASS.11) are at z f» 2.2 and 
the rest-frame is redshifted outside the Spitzer/MIPS photometric 

bands. An extrapolation would be uncertain and we choose not to estimate 
uLjj (60/im) for these sources. This plot is similar to that presented by Net- 
zer et al. (2007, their Fig. 5) and Lutz et al. (2008, their Fig. 5) and suggests 
that the star-formation rate of QSOs (measured by their 60/im luminosity) 
scales with AGN power, measured here by the 12/im luminosity. There is 
tentative evidence that 4/6 of the 2MASS QSOs plotted in the figure have 
more star-formation for their AGN power compared to optically selected 
PG QSOs, on average. 



sources were selected independent of their radio properties, unlike 
previous efforts to identify reddened QSOs, which started from ra- 
dio selected samples to identify objects with very red optical/near- 
IR colours (e.g. Gregg et al. 2002 ; White et al . 200 3 ; Glik man et all 
l2004l2007l : IUrrutiaet al. 2008a). 

Next we explore whether the radio detection rate of 2MASS 
QSOs is different from optically selected ones. For this exercise, 
we compare against the SDSS QSOs to take advantage of the large 
sample size and the well defined selectio n criteria. We use in par - 
ticular the SDSS-DR3 quasar catalog of I Schneider et al] ( l2005h , 
which was created by inspecting all spectra that were either tar- 
geted as quasar candidates, or classified as a quasar by the SDSS 
spectroscopic pipelines. This catalogue includes objects from all 
categories of spectroscopic target selection in the SDSS, not just 
those selected as QSO can didates. For our purpos es we use only 
targeted SDSS QSOs in the I Schneider et al] ( l2005h catalogue, i.e. 
sources that belong to the the main (or low-redshift) and the high- 
redshift QSO samples of the SDSS. These samples include QSOs 
with 15 < i < 19 mag that were sele cted based on the optical 
colour criteria of [Richards et al.l ( l2002h . We further apply a cut 
K < 14.5 mag to match the magnitude limit used to select 2MASS 
QSOs. There are 638 such sources in the Schneider et al. ( 200^ 
catalogue. The K = 14.5 mag magnitude cut does not introduce 
incompleteness in the sample, as it corresponds to i < 17.5 mag 
for the average (i — K) colour of optically selected QSOs, which 
takes values in the range 2-3, depending on redshift. Similarly, the 
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bright magnitude limit i > 15 mag for follow up spectroscopy of 
the SDSS QSO candidates corresponds to iC" > 12.5 mag. QSOs 
brighter than that may be missed from the SDSS sample. The sur- 
face density of such bright QSOs is likely to be small however, and 
it is not expected to introduce any significant bias in the SDSS op- 
tically selected QSO sample. 

In order to estimate the radio detection rate of the sample 
we use the FIRST radio surve y. There are 582 SDSS QSOs with 
K < 14.5 mag that fulfill the [Richards et alj C2002) criteria and 
overlap with the area covered by the FIRST survey. A total of 172 
of them, i.e. « 30 per cent, are detected at 1.4 GHz to the FIRST 
limit. This fraction should be compared with the detection rate of 
2MASS QSOs, i.e. 60 per cent. Using binomial statistics we esti- 
mate a 4 per cent probability of at least 6 radio detections in 10 trials 
given the success rate of 30 per cent of the optically selected QSO 
sample. The radio detection rate of 2MASS QSOs is higher than 
that of optically selected ones, albeit at the 96 per cent level (~ 2a). 
This is tentative evidence that red 2MASS QSOs are more often 
radio sources than optically selected QSOs. Glikman et alj ^20071) 
have identified reddened 2MASS QSOs by preselecting candidates 
in the FIRST radio survey. Their radio selected reddened sample 
represents a substantial fraction of the overall radio detected QSO 
population, up to 60 per cent. This finding further supports an as- 
sociation between red QSOs and radio emission. 

We caution that 2MASS QSOs are more luminous, on average, 
in the K-hsmd than the SDSS ones. This is demonstrated in Figure 
|7] plotting the Mk distribution of the two samples. Although the 
radio detection rate of SDSS QSOs does not appear to depend on 
Mk it might be possible that the somewhat different radio detec- 
tion rates of the two samples are because of 2MASS QSOs are more 
luminous. Larger samples of reddened QSOs are needed to further 
explore possible differences in their radio emission compared to 
optically selected samples. 



6.3 Fraction of reddened QSOs 

The fraction of reddened QSOs within the overall population is 

still under debate with different studies suggesting values from 15 

1^ II 1 

to over 5 per cent (e.g Richards et al. 2002; Wilkes et al. 2002; 

I White e t'al. 2003; Glikman et al. 2004). The difficulty in assessing 
this fraction is that one has to take into account the effect on dust in 
the detectability of reddened QSOs, which in effect means that at a 
given magnitude limit one detects only the most luminous reddened 
QSOs. In this section we take this effect into account by estimat- 
ing the maximum volume, Vmax, that a source could be detected 
taking into account the dust reddening and its luminosity. 

First however, we directly compare the surface density of 
the 2MASS QSOs with those of optically selected ones. In the 
previous section we cou n ted 63 8 QSOs in the SDSS-DR3 cata- 
log ue of [Schneider et al. I ( l2005l) with K < 14.5 mag that fulfill 
the [Richards et alj iiooj) criteria. Using the spectroscopic area of 
SDSS-DR3, 3732 deg^, we estimate a sky density of 0.17 deg"^ 
for optically selected QSOs lo K — 14.5 mag. In contrast, we find 
10 reddened 2MASS QSOs in the 5282 deg^ of the photometric 
SDSS-DR3 catalogue, which corresponds to a surface density of 
0.002 deg~^, 2dex lower than optically selected QSOs. The red- 
dened 2MASS QSOs selected here appear to be very rare and there- 
fore insignificant, in terms of numbers, for the overall population 
of QSOs. However, as shown in Figure|7]reddened QSOS are more 
luminous and the SDSS QSOs and it is therefore not surprising they 
are rarer. 

We address this issue using the 1/Vmax formalism (e.g. 
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Figure 7. distribution of optically selected SDSS QSOs with K < 
14.5 mag (open histogram) and reddened 2MASS QSOs (hatched his- 
togram). For the latter, is corrected for extinction. For both samples we 
plot only sources with redshift z < 1. This way we exclude high redshift 
lensed sources, such as 2MASS_11, with artificially enhanced luminosities. 
The inset plot shows as a function of Mj^ the fraction between the 2MASS 
and the total (2MASS -I- SDSS) QSO luminosity functions estimated us- 
ing the lA'max formalism as described in the text. At bright luminosities 
reddened 2MASS QSO represent 12% of the overall population. 



iLiUv et al.ll995h to estimate the if -band luminosity function of the 
2MASS QSOs relative to that of the overall QSO population. For 
each source with intrinsic absolute magnitude Mk the maximum 
volume Vmax is estimated given the selection function of the sam- 
ple, i.e. K < 14.5 mag, i? - if > 4 and J - AT > 2. For the calcu- 
lation of magnitudes, c olours and k-co rrections we adopt the SDSS 
composite QSO SED jVanden Berk et a l. 2001) and redden it by 
SMC type dust iPei 199^ assuming the optical extinction. Ay, 
listed in Table [3] We then sum the lA^max of each source in Mk 
luminosity bins. The same procedure is repeated for the optically 
selected SDSS QSOs with K < 14.5 mag, assuming Ay ~ for 
this population. In order to minimise evolution effects and to avoid 
high redshift sources for which lensing artificially boosts the ob- 
served luminosity (e.g. 2MASS_1 1), we only consider 2MASS and 
SDSS QSOs with z < 1. The results are shown in the inset plot of 
FigurelT] which plots as a function of AIk the ratio of the 2MASS 
sources luminosity function over the total luminosity function of 
both 2MASS and SDSS QSOs. The fraction of reddened QSOs is 
increasing from ~ 3 ± 2 per cent at Mk ~ —27.5 to about 12 ± 8 
per cent of the overall population at Mk = —29.5. This is evi- 
dence that reddened QSOs become important at bright luminosities. 
We caution that small number statistics are an issue in this calcula- 
tion. Larger reddened QSO samples at magnitudes fainter than the 
2MASS K = 14.5 limit are needed to further explore this trend and 
to provide a more reliable measure of the reddened QSO fraction. 
Efforts to identify QSOs in the the UKIRT Infrared Deep Sky Sur- 
vey (UKIDSS, K ~ 18 mag) incl uding dust reddened ones, have 
recently appeared in the literature jMaddox et al.ll2008h . although 
the source selection is very different to that presented here. 
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7 DISCUSSION 

This paper combines photometric data from the largest optical 
and near-IR surveys, SDSS and 2MASS, to select QSOs with red 
optical/near-IR colours, which are atypical to those of the tradi- 
tional population of optically selected broad-line QSOs. It is shown 
that the UV/optical SEDs of the red QSOs can be explained by 
moderate amounts of dust extinction, Av = 1.3 — 3.2. An open 
question is whether red QSOs can be explained by orientation ar- 
guments, and therefore are a subset of optically selected QSOs, or 
they represent a different stage of the life of QSOs. 

On the latter point, there is a large volume of literature on the 
possible evolution link between starbu rsts, Q SOs and spheroids 
Both observations (e . g. Sanders et^ 1 19881: iSanders & Mirabel 
1994 IClements|[2OO0l: ICanalizo & Stocktonll200ll; iKomossa et al 



2003 : [Alexander et al."2005l: IVeilleux et al.ll2006l : ISchweitzer et al 
20061 : iDasvraetal. 20oJ_ and numerical simulations of ma- 



jor mergers (e.g. JPiMatteoetalJ l2005l : ISpringel et al] l2005l : 
IHopkins et al.ll20()3) support such an association. It is proposed in 
particular, that some mechanism, possibly mergers, drive gas effi- 
ciently to the central galaxy regions, thereby triggering powerful 
starbursts and black hole growth at a high accretion rate. During 
this stage both the central engine the star-forming regions are en- 
shrouded in dust and gas clouds. As a result the system appears 
luminous in the infrared (LIRGs or ULIRG stage). AGN and/or 
starburst driven winds (feedback) will eventually develop and blow 
away the dust cocoon at the nuclear regions, thereby allowing the 
AGN to shine unobscured. This is the stage where the system can 
be identified as UV/optically luminous broad line QSO. When all 
the cold gas in the nuclear regions is depleted, either because it 
is used up in star-formation, consumed by the black-hole, heated 
or pushed away by the feedback processes, the AGN and starburst 
will switch off, leaving behind a quiescent galaxy. In this picture 
reddened QSOs are systems captured just before or during the final 
blow-out stage of their evolution. 

Our findings, although limited by the small number of sources 
in the sample, are in broad agreement with this scenario. The far-IR 
properties of about 70 per cent of the reddened QSOs are consistent 
with starburst activity. There is also evidence that the star-formation 
rate of these systems, normalised to the AGN power, is higher than 
that of optically selected QSOs, on average. A possible interpre- 
tation of this finding is that feedback from the central engine has 
not yet suppressed the star-formation in the host galaxy at the same 
level as in typical QSOs. There is also tentative evidence, signif- 
icant at the 96 per cent confidence level, that reddened QSOs are 
more often associated with radio sources compared to optically se- 
lected QSOs. If this is confirmed by larger samples, it may suggest 
the onset of the AGN driven feedback in these systems via radio 
jets. The tentative evidence for larger dust covering fractions (ap- 
proximated by Ltor/Lopt) of the 2MASS QSOs presented in this 
paper compared to optically selected SWIRE QSOs, is also con- 
sistent with the young AGN interpretation. With the data presented 
here we cannot com ment on the o pt ical m orphology of reddened 
QSO host galaxies. iHutchings et alj ( l2003l.l200^ and lUrrutia et al.l 
t2008b,) explored this issue using a different sample of red 2MASS 
QSOs selected to have near-IR colours J—K > 2 mag. They report 
a higher incidence of systems with interaction signs in their sampl e 
compared to optically selected QSOs (but see lMarble et alj|2003l) . 
Although this conclusions is limited by the absence of an appro- 
priate comparison sample, it is consistent with the merger induced 
starburst/AGN activity picture above, and suggests that reddened 
QSOs are at an earlier evolutionary stage than optically selected 



ones. The relatively low fraction of 2MASS QSOs, i.e. 3-12 per 
cent of the overall QSO population depending on Mk, may there- 
fore reflect the short timescale of the particular evolutionary stage 
relative to the QSO lifetime. 

Another class of sources that are also proposed to pre- 
cede the optically luminous QSO p hase are ULIRGs with broa d 
emission lines (type-I ULIRGs: e.g. iKawakatu et al.ll2003l. l2006h . 
These systems have young stellar populations (< 300 Myr old 
ICanalizo & StocktonI 2001), a large fraction of them experienc e 
tidal interactions jCanalizo & Stocktoij|200ll : iLipari et alJboOSh , 
their black hole masses are smaller compared to ellipticals and op- 
tically selected QSOs for the size of their bulge ( Kawaka tu et al.l 
|2006), Broad Absorption L ines (BAL) indicative of o utflows are 
frequent in this population jCanalizo & Stocktonll2001i) , and most 
of them have narrow permitted emission lines (e.g. FWHM of 
H/3 less than 4000 km s~^) similar to narrow-line Seyfert 1 galax- 
ies, which are also believed to be young AGN (Zheng et al. 20o3: 
lAnabukilliooi) . The 2MASS QSOs studied here share some of 
the properties of type-I ULIRGs, i.e. Lm > 10^^ LQ, permit- 
ted emission lines with t<WHM typically < 2000 km s~^, ongo- 
ing star-formation. Also, IU rrutiaetal] ( l2008 i^ have recently found 
an unusually high fraction of BAL sources among red 2MASS 
QSOs selected with criteria similar to those used here. It is there- 
fore tempting to link the two types of sources. A major differ- 
ence between the two populations however, is that Type-I ULIRGs 
are typically less reddened at UV/optical than 2MASS QSOs (e.g. 
Canalizo & Stockton 2001). Is it therefore possible that 2MASS 
QSOs are even younger than type-I ULIRGs in the evolution 
scheme outlined above? 

Figure [8] explore this possibility by comparing the infrared 
colours, a(60, 25) against q(100, 60) (see figure caption for de- 
tails on the estimation of the colours), of PG and 2MASS QSOs 
with the sample of type-I ULIRGs of Canalizo & Stockton (2001). 
Different types of sources occupy fairly distinct regions in this fig- 
ure. Sources dominated by cool dust are found in the left part of 
the diagram, close to line marked "black body". Systems domi- 
nated by non-thermal emission populate the left part of the dia- 
gram, close to the "power-law" line. According to the evolution 
scenario described above QSOs should evolve from right to left 
in Figure [S] Optically selected PG QSOs dominate the upper left 
part of the diagram, the type-I ULIRGS of Canalizo & Stockton 
(2001) are selected to lie between the two diagonal lines, while the 
four low redshift 2MASS QSOs (z < 0.6) scatter toward the far 
right end of the diagram. This suggests, but does not prove, that 
2MASS QSOs may represent young QSOs, possibly before they 
become type-I ULIRGs and eventually typical UV/optically lumi- 
nous broad line QSOs. Clearly, more data are needed to further 
explore this scenario. What fraction of 2MASS QSOs qualify for 
narrow line Seyfert Is? What are the black hole masses and the ac- 
cretion rates of these systems? Do we see optical spectral evidence 
for outflows (e.g. abso rption features), similar to those observed 
in type-I ULIRGs (e.g. lCanalizo & Stocktorj|200lh ? What are the 
properties of their host galaxies? All these are questions that need 
to be addressed before concluding on the nature of 2MASS QSOs. 



8 CONCLUSIONS 

We combine optical and near-IR photometry from the ground with 
Spitzer mid- and far-IR broadband data to explore the SEDs of a 
sample of 10 QSOs with red optical and near-IR colours selected 
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a(100,60) 

Figure 8. Infrai'ed colour-colour plot. The colour indices ai'e defined 
a(Ai,A2) = -log(/Ai//Ai)/log{Ai/A2), where fx is the flux den- 
sity at wavelength A. The filled circles are the sample of PG QSOs used 
here, open circles are the type-I ULIRGS of Canalizo & Stockton (2001) 
and the (red) crosses are the low redshift 2MASS QSOs (z < 0.6) in the 
present sample (total of 4). For higher redshift 2MASS QSOs at least one of 
the rest-frame wavelengths 25, 60 or 100/.tm, is outside the Spitzer spectral 
window. For the 2MASS and PG QSOs the rest-frame flux densities used to 
estimate the infrared colour indices are estimated by interpolating between 
the data points of the observed SEDs. The infrared colours of the Canal- 
izo & Stockton (2001) sample are estimated using the IRAS flux densities. 
The diagonal li nes show the position of power-law and black-body SEDs 
(adapted from iLiparill994l) . 



by combining the 2MASS with the SDSS. The main conclusions 
are 

• The UV/optical part of the SED can be fit assuming an in- 
trinsic spectrum similar to that of typical optically selected QSOs 
reddened by moderate amounts of dust. 

• There is evidence that red 2MASS QSOs have higher level 
of star-formation (approximated by the 60^m luminosity) for their 
AGN power (measured by the 12/im luminosity), compared to op- 
tically selected PG QSOs. 

• We find that 2MASS QSOs are more often associated with 
radio emission, compared to optically selected SDSS QSOs. This 
excess, although interesting, is only significant at the 96% level 

• The fraction of red QSOs in the overall population increases 
from 3 ± 2% at Mk = -27.5 mag to 12 ± 8% at Mk = 
—29.5 mag, suggesting that these systems become more important 
at the bright luminosities. 

• We argue that these findings are consistent with a pic- 
ture where reddened QSOs are young AGN shortly before they 
shed their cocoon of dust and gas clouds and become typ- 
ical UV/optically luminous QSOs. Comparison with Seyfert-1 
ULIRGs, another class of sources proposed to be young QSOs, sug- 
gests that red 2MASS QSOs likely represent an earlier evolutionary 
stage. 
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APPENDIX A: NOTES ON INDIVIDUAL SOURCES 

2MASS QSO#2: The optical spectrum of this sources shows broad 
Balmer emission lines, H/3, H7. Fitting a Gaussian to the H/3 we 
estimate a FWHM of ^ 2000kms"\ The [Nelll] 3869 line is 
also seen in the spectrum with broad profile. It is interesting that 
the Balmer lines are redshifted from the [OIII] 4959 and 5007A 
lines by ~ 1500 kms~^. A similar effect was reported by Mcin- 
tosh et al. (1999) for luminous high redshift QSOs. In their sam- 
ple the mean redward offset of the H beta line with respect to the 
[OIII] 5007 was estimated to be ~ 520km s^^. Visual inspection 
of the optical spectrum also suggests that the Balmer lines have 
asymmetric profiles with blue wings. 

2MASS QSO#3: This source has narrow emission lines, H/3, 
[OIII] 4959 and 5007A [NII]6583A, and the Ha line with broad 
wings. A simultaneous two component Gaussian fit to the blended 
Ha and [Nil] 6583A lines fails to reproduce the flux in the wings 
of Ha. An additional gaussian is introduced to account for the 
broad wings. The FHWM of this component is estimated to be 
~ 1500 kms~^. We also estimate diagnostic emission line ratios 
using the H/3, [OIII] 5007A [Nil] 6583A, and the narrow com- 
ponent of the Ha. We estimate log([OIII]/H/3) = +0.50 and 
log([NII]/Ha) = —0.58, which place this s ource in the Transitio n 
Objects region of the diagnostic diagram of Kewley et al. ( 200lj). 
The radio emission from th is source is con sistent with the starburst 
radio/far-IR correlation of ICondoi] jl992l) . Adopting the SED of 
M 82 to extrapolate the measured rest-frame 60 fim flux density to 
100 /im we estimate q — 2.19, which is consistent with the media n 
for starbursts g = 2.3 ± 0.2 ( iHelou et alJI 19851 : ICondorilll992h . 
We estimate an extinction of Ay = 1.6 from the Balmer decre- 
ment of the narrow emission line components. This compares well 
to Av = 1.9 determined from the SED fits to the continuum. 
2MASS QSO#4: Only narrow emission lines are observed in the 
optical spectrum of this source, H/3, [OIII] 4959 and 5007A. The 
spectral window does not cover the Ha or the Paschen line region 
which suffer less extinction. This source is therefore likely to also 
have broad lines. 

2MASS QSO#5: This source has optical spectroscopic observa- 
tions, which however did not yield a redshift determination. SED 
fitting estimates a photometric redshift z = 3.1. 
2MASS QSO#6: We find broad Ha, Paa and Pa/3 emission lines. 
A three component Gaussian fit to the Ha line (to account for the 
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blended [Nil] 6583 A line) gives a FWHM of 2200 kms^^ for 
the broad component. The FWHM of the Paschen lines is about 
1500 kms~^, estimated by single Gaussian fits to the line profiles. 
The extinction is Av = 5.7 from the Balmer decrement. In this cal- 
culation we use the narrow line Ha component to compare with the 
narrow H/9 line. The extinction above is higher than Av ~ 2.3 de- 
termined from the SED fits to the continuum. The Paschen line ratio 
is estimated Paa/Pa/S = 1.9 and is consistent with Av = 5.3 or 
Av = 9.6 adopt ing respectively, a n empiri cally determined intri n- 
sic ratio of 1.05 dSoifer et al.ll2004h or 0.64 jOlikman et alj2006l) . 
2MASS QSO#7 : Optical spectroscopic observations are available, 
but did not yield a redshift determination. SED fitting estimates a 
photometric redshift z = 2.73. 

2MASS QSO#8: This source lies &t z = 2.16 and has broad Ha 
with FWHM of about 4700 kms~^, estimated by single Gaussian 
fit to the line profile. The extinction is Av = 4.4 from the Balmer 
decrement. In this calculation we use single guassian line fits to the 
Ha and H/3 lines. The extinction above is higher than the Av ~ 
1.2 determined from the SED fits to the continuum. 
2MASS QSO#10: The optical spectrum of this source shows 
narrow emission lines only, [Oil] 3727A, H/3, [OIII]4959 and 
5007Athe near-IR spectrum of this source has been presented by 
Glikman et al. (2004) and shows broad Pa/3 emission line. 
2MASS QSO#ll: This is the lensed reddened QSO first found by 
Gregg et al. (2002). The optical and near-IR spectrum presented in 
that paper shows broad emission lines. 

2MASS QSO#12: This source has broad Paa and Pa/3 emis- 
sion lines with FWHM of about 1600 kms^^, estimated by single 
Gaussian fits to the line profiles. The Paschen line ratio is estimated 
Paa/Pa/3 — 1.7. The corresponding reddening is Av ~ 4.3 
or A v = 8.6 for the e mpirically determined intrinsic ratios of 
1.05 jSoiferetalJl2004h and 0.64 CGlikman et al.. .2006) respec- 
tively. These ratios are larger than the reddening of Av ~ 3 de- 
termined from the SED fits to the continuum. This source has also 
been selected in the sample of Glikman et al. (2007) and Urrutia et 
al. (2008a). 
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